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Chapter 1 
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Alex PomaroQ 

Departament de Fisica, Universitat Autbnoma de Barcelona, 
08193 Bellaterra, Barcelona 

Andrea WulzeiQ 

Institut de Theorie des Phenomenes Physiques, EPFL, 
CH-1015 Lausanne, Switzerland 

We review the procedure to calculate baryonic properties using a recently pro- 
posed five-dimensional approach to QCD. We show that this method give pre- 
dictions to baryon observables that agree reasonable well with the experimental 
data. 

1.1. Introduction 

In 1973 Gerard 'tHooft proposed, in a seminal article [1], a dual description for QCD. 
He showed that in the limit of large number of colors (N c ) strongly-interacting gauge 
theories could be described in terms of a weakly-interacting theory of mesons. It 
was later recognized [2] that, in this dual description, baryons appeared as solitons 
made of meson fields, as Skyrme had pointed out long before [3]. These solitonic 
states were therefore referred as skyrmions. 

Skyrmions have been widely studied in the literature, with some phenomenolog- 
ical successes [4]. Nevertheless, since the full theory of QCD mesons is not known, 
these studies have been carried out in truncated low-energy models either incorpo- 
rating only pions [2, 3] or few resonances [4]. It is unclear whether these approaches 
capture the physics needed to fully describe the baryons, since the stabilization 
of the baryon size is very sensitive to resonances around the GeV. In the original 
Skyrme model with only pions, for instance, the inverse skyrmion size p~ l equals 
the chiral perturbation theory cut-off A x pt ~ ^itF^ (as it should be, since this is the 
only scale of the model), rendering baryon physics completely incalculable. Other 
examples are models with the p-meson [5] or the w-meson [6] which were shown 
to have a stable skyrmion solution. The inverse size, also in this case, is of order 
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Trip ~ A x px, which is clearly not far from the mass of the next resonances. Includ- 
ing the latter could affect strongly the physics of the skyrmion, or even destabilize 
it. 

In the last ten years the string/gauge duality [7-9] has allowed us to gain new 
insights into the problem of strongly-coupled gauge theories. This duality has been 
able to relate certain strongly-coupled gauge theories with string theories living in 
more than four dimensions. A crucial ingredient in these realizations is a (compact) 
warped extra dimension that plays the role of the energy scale in the strongly- 
coupled 4D theory. This has suggested that the QCD dual theory of mesons pro- 
posed by 'tHooft [1] must be a theory formulated in more than 4 dimensions. 

Inspired but this duality, a five-dimensional field theory has been proposed in 
Refs. [10, 11] to describe the properties of mesons in QCD. This 5D theory has a 
cut-off scale A5 which is above the lowest-resonance mass m p . The gap among these 
two scales, which ensures calculability in the meson sector, is related to the number 
of colors N c of QCD. In the large iV c -limit, one has A^/nip — ► 00 and the 5D model 
describes a theory of infinite mesonic resonances, corresponding to the Kaluza-Klcin 
(KK) spectrum. This 5D model has provided a quite accurate description of meson 
physics in terms of a very limited number of parameters. 

Further studies, boosted by this success, have recently shown that the 5D model 
can also successfully describe baryon physics [12-14]. As Skyrme proposed [3], 
baryons must appear in this 5D theory as solitons. These 5D skyrmion-likc solitons 
have been numerically obtained (see Fig. ll.l|) and their properties have been studied. 
Their inverse size pj 1 ~ m p have been found to be smaller than the cut-off scale 
A5, showing then that, contrary to the 4D case, they can be consistently studied in 
5D effective theories. Indeed, the expansion parameter which ensures calculability 
is provided by l/(p s A 5 ) -C 1. 




Fig. 1.1. Energy density, in the plane of the 4D radial and the extra fifth coordinate, of the 
skyrmion in a 5D model for QCD. 



In this article we will review the properties of baryons obtained in Refs. [12-14] 
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using the five-dimensional model of QCD of Rcfs. [10, 11, 15]. We will show how 
the calculation of the static properties of the nucleons, such as masses, radii and 
form factors, are performed, and will compare the predictions of the model with 
experiments. As we will see, these predictions show a reasonably good agreement 
with the data. 

There have been alternatives studies to baryon physics using 5D models. Nev- 
ertheless, these studies have encountered several problems. For example, the first 
approaches [16] truncated the 5D theory and only considered the effects of the first 
resonances. This leads to skyrmions whose size is of the order of the inverse of the 
truncation scale, and therefore sensitive to the discarded heavier resonances. Later 
studies [17-19] were performed within the Sakai-Sugimoto model [20]. It was shown, 
however, that baryons are not calculable in this framework as their inverse size is 
of the order of the string scale which corresponds to the cut-off of the theory [17]. 

1.2. A five-dimensional model for QCD mesons 

The 5D model that we will consider to describe mesons in two massless flavor 
QCD is the following. This is a U(2)l x U(2)r gauge theory with metric ds 2 = 
a(z) 2 {j]^ u dx^dx v — dz 2 ), where x 11 represent the usual 4 coordinates and z, which 
runs in the interval [zuv, z m)i denotes the extra dimension. We will work in AdS5 
where the warp factor a(z) is 

= (1-1) 

z 

and ztjv ^ to be taken at the end of the calculations. In this limit zi R coincides 
with the AdS curvature and the conformal length 

L= dz. (1.2) 

J ZTJV 

The U(2)l and U(2)r gauge connections, denoted respectively by Lm and Rm 
(M = {/z, 5}), are parametrized by Lm = L M <r /2 + and Rm = R < MO' a /2 + 

-RmI/2, where a a are the Pauli matrices. This chiral gauge symmetry is broken by 
the conditions on the boundary at z = Zir (IR-boundary) , which read 

(L p - R M ) |z=z m = , (L M 5 + R^s) \ z =z m = , (1.3) 

where the 5D field strength is defined as Tjmn = 9m^n — c^vLm — i\Lm, Ljv]) 
and analogously for Rmjv- At the other boundary, the UV one, we can consider 
generalized Dirichlet conditions for all the fields: 

L^x |z=z uv = Ifi j R^j |z=zuv = r ti ■ (1-4) 

The 4D fields Z M and r M are arbitrary but fixed and they can be interpreted, as we 
will now discuss, as external sources for the QCD global currents. We will eventually 
be interested in taking the sources to vanish. 
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We can now, inspired by the "holographic" formulation of the AdS/CFT corre- 
spondence [7-9], try to interpret the above 5D model in terms of a 4D QCD-like 
theory, whose fields we will generically denote by \&(x) and its action by S4. This 
is a strongly coupled 4D theory that possesses an U(2)l x U(2)r global symmetry 
with associated Noether currents R . If the 4D theory were precisely massless 
QCD with two flavors, the currents would be given by the usual quark bilinear, 

\3l r) = Q'l rI^Ql r- Defining Z[l u , r u ] as the generating functional of current 

\ ' / ij 

correlators, we state our correspondence as 



X»J>cxp 



J VL M VR M exp [iS B [L, R]] 



(1.5) 



where the 5D partition function depends on the sources lp, r M through the UV- 
boundary conditions in Eq. (| 1 .4[) . 

Eq. (|1.5[) leads to the following implication. Under local chiral transformations, 
Z receives a contribution from the U(2) 3 anomaly, which is known in QCD. 0This 
implies [9, 21, 22] that the 5D action must contain a Chern-Simons (CS) term 



S, 



cs 



24tt 2 



Au; 5 (L)-a; 5 (R)] 



(1.6) 



whose variation under 5D local transformations which does not reduce to the iden- 
tity at the UV exactly reproduces the anomaly. The CS coefficient will be fixed to 
N c = 3 when matching QCD. The CS 5-form, defining A = —iAudx , is 



w 5 (A) = Tr 



A(dA) 2 + -A 3 (dA) + -A 5 
2 5 



(1.7) 



When A is the connection of an U(2) group, as in our case, one can use the fact 
that SU(2) is an anomaly-free group to write W5 as 



w 6 (A) = ^ATr [F 2 ] + ±A (dA 



dTr 



AAF - -AA 3 

4 



(1.8) 



where A = A + Al/2 and A is the SU(2) connection. The total derivative part 
of the above equation can be dropped, since it only adds to Scs an UV-boundary 
term for the sources. 

The full 5D action will be given by S5 = S g + Scs , where S g is made of locally 
gauge invariant terms. S g is also invariant under transformations which do not 
reduce to the identity at the UV-boundary, and for this reason it does not contribute 



a The 5D semiclassical expansion we perform in our model corresponds, as we will explain in the 
following, to the large- N c expansion. This is why we are ignoring the U(1)-SU(N C ) 2 anomaly of 
QCD, which is subleading at large- N c . Being this anomaly responsible for the if mass, our model 
will contain a massless rj' . 
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to the anomalous variation of the partition function. Taking the operators of the 
lowest dimensionality, we have 

S B = -J d 4 xj* m dza(z) jlr [L M nL mn ] + ^L MN L MN + {L «-> R} 

(1.9) 

We have imposed on the 5D theory invariance under the combined {x — > — x, L <-> 
R}, where x denotes ordinary 3-space coordinates. This symmetry, under which 
Scs is also invariant, corresponds to the usual parity on the 4D side. We have 
normalized differently the kinetic term of the SU(2) and U(l) gauge bosons, since 
we do not have any symmetry reason to put them equal. In the large- N c limit of 
QCD, however, the Zweig's rule leads to equal couplings (and masses) for the p 
and w mesons, implying a = 1 in our 5D model. Since this well-known feature of 
large- N c QCD does not arise automatically in our 5D framework (as, for instance, 
the equality of the p and ui masses does), we will keep a as a free parameter. The 
CS term, written in component notation, will be given by 

Scs = J^J d 5 xlle M » op QL M Tr[L N oL PQ ] 

+ ^ MNOPQ L M L NO L PQ - {L ~ R}\ . (1.10) 

The 5D theory defined above has only 3 independent parameters: M5, L and a. 

Let us make again use of Eq. (|1.5[) to determine the current operators through 
which the theory couples to the external EW bosons. These currents are obtained 
by varying Eq. (|1.5p with respect to Z M (exactly the same would be true for r M ) and 
then taking Z M = = 0. The variation of the l.h.s. of Eq. (|1.5[) simply gives the 
current correlator of the 4D theory, while in the r.h.s. this corresponds to a variation 
of the UV-boundary conditions. The effect of this latter can be calculated in the 
following way. We perform a field redefinition L M — ► L M + <5L M where Sh fl (x 1 z) is 
chosen to respect the IR-boundary conditions and fulfill Sli^(x, zuv) = 51^. This 
redefinition removes the original variation of the UV-boundary conditions, but leads 
a new term in the 5D action, 8S§. One then has 

i J d 4 xTr [(j1(x))SI^x)} = 1 J VL M VTL M 5S 5 [L, R] exp [tS 5 [L, R]] , (1.11) 

where the 5D path integral is now performed by taking = = 0, i.e. normal 
Dirichlet conditions. The explicit value of 5S§ is given by 

8S 5 = J d 4 xTT[Jl(x)Sl tl (x)} + J d 5 s(EOM) • SL , (1.12) 

where 3 L ^ = J c l^a a + J Lll t and 

Jl, = M 5 (a(z)L; 5 ) \ z=zvv , J iM = a 2 M 5 (a(z)L^ 5 ) | 2=Zuv . (1.13) 

The last term of Eq. (|1.12j) corresponds to the 5D "bulk" part of the variation, which 
leads to the equations of motion (EOM) . Remembering that the EOM always have 
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zero expectation value El, we find that we can identify J£ of Eq. (|1.13|) with the 
current operator on the 5D side: O'^Kd = (J^sd- Notice that the CS term has 
not contributed to Eq. (|1.12|) due to the fact that each term in Scs which contains 
a d z derivative (and therefore could lead to a UV-boundary term) also contains L M 
or R M fields; these fields on the UV-boundary are the sources and ry, that must 
be put to zero. 

1.2.1. Meson Physics and Calculability 

The phenomenological implications for the lightest mesons of 5D models like the 
one described above have been extensively studied in the literature. Let us briefly 
summarize the main results here. If rewritten in 4D terms, the theory contains 
massless Goldstone bosons that parametrize the U(2)l x U(2)r/U(2)v coset and 
describe the pion triplet and a massless rf . The pion decay constant is given by 

The massive spectrum consists of infinite towers of vector and axial-vector spin- 
one KK resonances. Among the vectors we have an isospin triplet, the p^ n \ and a 
singlet uj( n \ The axial- vectors are again a triplet and a singlet f[ n K We want 
to interpret, as our terminology already suggests, the lightest states of each tower 
as the jo(770), w(782), ai(1260) and /i(1285) resonances, respectively. The model 
predicts at leading order, i.e. at tree- level, 

m p = m u ~ — , m ai =m fl ~—, (1.15) 

compatibly with observations. The model also predicts the decay constants Fi 
and couplings gi for the mesons as a function of Afs, L and a that can be found 
in Refs. [10-12, 15]; here we only notice, for later use, their scaling with the 5D 
coupling: 

F~y/W 5 , 3*-^=, (1-16) 

V M 5 

while the masses, as shown above, do not depend on M5. In Table ITTTI we show a fit 
to 14 meson quantities. The best fit is obtained for the values of 1/L = 343 MeV, 
M5I/ = 0.0165 and a = 0.94 for the three parameters of our model. The minimum 
Root Mean Square Error (RMSE) corresponding to those values is found to be 11% 
and the relative deviation of each single prediction is below around 15%. 

Concerning the choice of the meson observables, some remarks are in order. 
First of all, we are only considering the lowest state of each KK tower because we 
expect the masses and couplings of the heavier mesons to receive large quantum 
corrections. Our model is indeed, as we will explain below, an effective theory valid 

b We have actually shown this here; notice that <5L,j was completely arbitrary in the bulk, but the 
variation of the functional integral can only depend on <5/ M = <5L fl (x, ztjv)- 
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Table 1.1. Global fit to mesonic physical quantities. 
Masses, decay constants and widths are given in MeV. 
Physical masses have been used in the kinematic factors 
of the partial decay widths. 

Experiment AdSs Deviation 



m p 


775 


824 


+6% 


m ai 


1230 


1347 


+10% 


m u 


782 


824 


+5% 


F P 


153 


169 


+11% 


F^/F p 


0.88 


0.94 


+7% 


Ftt 


87 


88 


+1% 




6.0 


5.4 


-10% 


Lg 


6.9 ■ 10~ 3 


6.2 ■ 10~ 3 


-10% 


Lio 


-5.2 ■ 10" 3 


-6.2 ■ 10" 3 


-12% 


T(uj — » 7T7) 


0.75 


0.81 


+8% 


r(w -> 3tt) 


7.5 


6.7 


-11% 


r(p -> 7T7) 


0.068 


0.077 


+13% 


r(o; — ► TTpp) 


8.2 ■ 10~ 4 


7.3 ■ 10- 4 


-10% 


F(lu — » 7ree) 


6.5 ■ 10~ 3 


7.3 ■ 10~ 3 


+12% 



up to a cut-off A5 ~ 2 GeV and our tree-fevel calculations only correspond to the 
leading term of an E/A5 expansion. Apart from this restriction, we must include in 
our fit observables with an experimental accuracy better than 10%. This is because 
we want to neglect the experimental error in order to obtain an estimate of the 
accuracy of our theoretical predictions. Much more observables can be computed, 
once the best-fit value of the parameters are obtained, and several of them have 
already been considered in the literature. For instance, one can study the other 
low-energy constants of the chiral lagrangian, the physics of the /1 resonance or 
the pseudo-scalar resonances which arise when the explicit breaking of the chiral 
symmetry is taken into account [11]. It would also be interesting to compute the 
ai — > 7T7 decay, which is absent in our model at tree-level and only proceeds via 
loop effects or higher-dimensional terms of our 5D effective lagrangian. El 

As discussed in the Introduction, the semiclassical expansion in the 5D model 
should correspond to the large- N c expansion on the 4D side. The results presented 
above provide a confirmation of this interpretation: at large- N c meson masses are 
expected to scale like AT?, while meson couplings and decay constants scale like 
gi,l/Fi ~ l/yQVc- These scalings agree with Eq. (|1.15[) and () 1 . 1 6(1 if the parameters 



igrei 

a, L and M 5 are taken to scale like|£| 

a~N° , L~N° , M 5 ~JV C . (1.17) 
This leads us to define the adimcnsional iV c -invariant parameter 

7^ 1C J!° , , (1.18) 



c Higher order contributions will also change our tree-level prediction Lg + Lio 
related with the absence of the ai-n-'y vertex. 

d This scaling can also be obtained from the AdS/CFT correspondence. 



= 0, which is again 
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whose experimental value is 7 = 1.23 and will be useful later on. We will also show 
in the following that the assumed scaling of the 5D parameters leads to the correct 
N c scaling in the baryon sector as well. 

Other descriptions of vector mesons in terms of massive vector fields, i.e. models 
with Hidden Local Symmetry (HLS) [23, 24] or two-form fields [25], also correctly 
reproduce the meson physical properties. Nevertheless,we believe that 5D models, 
as the one discussed here, present more advantages. □ First of all, they contain 
less parameters. In the models of Refs. [23-25], for example, the mass and the 
couplings of each meson are independent parameters; also anomalous processes, 
those involving an odd number of pions, depend on several operators with unknown 
coefficients which arise at the same order, while in our case they all arise from a 
single operator, the 5D CS term. Finally, Vector Meson Dominance is automatic in 
our scenario, while it needs to be imposed "by hand" in the case of HLS. 

Moreover, and perhaps more importantly, 5D models are calculable effective 
field theory in which higher-dimensional operators are suppressed by the cut-off of 
the theory A5. Calculations can be organized as an expansion in E/A5, where E 
is the typical scale of the process under consideration. Given that the cut-off is 
parametrically bigger than the mass of the lightest mesons, reliable calculations of 
masses and couplings can be performed. 

Let us now use naive dimensional arguments to estimate the maximal value of 
our cut-off A5. This is determined by the scale at which loops are of order of tree- 
level effects. Computing loop corrections to the F 2 operator, which arise from the 
F 2 term itself, one gets A5 ~ 247r 3 Ms. Nevertheless, one gets a lower value for 
A 5 from the CS term due to the N c dependence of its coefficient. Indeed, at the 
one-loop level, the CS term gives a contribution of order M5 to the F 2 operator for 
A5 ~ 24ir 3 M§/Nc ^ . Even though the cut-off scale lowers due to the presence of 
the CS term, we can still have, in the large- N c limit, a 5D weakly coupled theory 
where higher-dimensional operators are suppressed. The cut-off can be rewritten as 

A 5 ~^^V~2GeV, 
2 7«L 

where we have used the best-fit value of our parameters. 

The power of calculability of our 5D model makes it very suitable for studying 
baryon physics. Indeed, the typical size of the 5D skyrmion solution will be of order 
p s ~ 1/rrip, guaranteeing that effects from higher-dimensional operators will be 
suppressed by m p /A 5 ~ 0.4. This is therefore, we believe, the first fully consistent 
approach towards baryon physics. 



e It must be possible, generalizing what was done in Ref. [26], to rewrite our model as a 4D HLS 
with infinitely many U (2) hidden symmetry groups. The comparison with HLS models that we 
perform in this section only applies, therefore, to the standard case of a finite number of hidden 
symmetries. 
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1.3. Baryons from 5D Skyrmions 
1.3.1. 4D Skyrmions from 5D Solitons 

Time-independent configurations of our 5D fields, which correspond to allowed ini- 
tial (t — > — oo) and final (f — * +00) states of the time evolution, are labeled by the 
topological charge 

1 f rzm 

B = — d 3 x dze^Tr [L^L^ - BTRi"] , (1.19) 

J Jz vv 

where the indeces (i,0, . . . run over the 4 spatial coordinates, but they are raised 
with Euclidean metric. We will now show that B can only assume integer values, 
which ensures that it cannot be changed by the time evolution. This makes B a 
topologically conserved charge which we identify with the baryon number. In order 
to show this, and with the aim of making the relation with the skyrmion more 
precise, it is convenient to go to the axial gauge L 5 = R 5 =0. The latter can 
be easily reached, starting from a generic gauge field configuration, by means of a 
Wilson-line transformation. In the axial gauge both boundary conditions Eqs. (|1.3p 
and (|1.4|) (in which we take now I = r = 0) cannot be simultaneously satisfied. Let 
us then keep Eq. (|1.3|) but modify the UV-boundary condition to 



U U= zuv = i tf(x)fi(l7(x)t , R t \ z=zuv = , (1.20) 

where Li and Ri are the gauge fields in the axial gauge and i runs over the 3 ordinary 
space coordinates. The field U (x) in the equation above precisely corresponds to 
the Goldstone field in the 4D interpretation of the model [22] . Remembering that 
F A F = du>3, where is the third CS form, the 4D integral in Eq. (|1.19jl can be 
rewritten as an integral on the 3D boundary of the space: 

B = ^2 / [«s(Z)- "3(5)1 . (1.21) 
8tt 2 J 3D L J 

The contribution to B coming from the IR-boundary vanishes as the L and R terms 
in Eq. (|1.21[) cancel each other due to Eq. (|1.3|l . This is crucial for B to be quantized 
and it is the reason why we have to choose the relative minus sign among the L 
and R instanton charges in the definition of B. At the x 2 — ► 00 boundary, the 
contribution to B also vanishes since in the axial gauge d$Ai = (in order to have 
Fbi — 0). We are then left with the UV-boundary which we can topologically regard 
as the 3-sphere S3. Therefore, we find 



uv 



L l {=iUd l U ] ) 



1 



24tt 2 



d 3 xe ijk Tr [Uditf UdjU* Udktf] GZ. (1.22) 



The charge B is equal to the Cartan-Maurer integral invariant for SU(2) which is 
an integer. 
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In the next section we will discuss regular static solutions with nonzero B. If 
they exist, they cannot trivially correspond to a pure gauge configuration. Moreover, 
the particles associated to solitons with B = ±1 will be stable given that they have 
minimal charge. Eq. (|1.22p also makes the relation with 4D skyrmions explicit: 
topologically non-trivial 5D configurations are those for which the corresponding 
pion matrix U (x) is also non-trivial. The latter corresponds to a 4D skyrmion with 
baryon number B. In a general gauge, the skyrimion configuration J7(x) will be 
given by 

t/(x) = P jcxp -i [ dz'R 5 (x,z') j-p[exp i f dz'L 5 (x,z') X , (1.23) 

where P indicates path ordering. From a 4D perspective, the 5D soliton that we 
are looking for can be considered to be a 4D skyrmion made of Goldstone bosons 
and the massive tower of KK gauge bosons. 



1.3.2. The Static Solution 

In order to obtain the static soliton solution of the 5D EOM of our theory it is 
crucial to specify a correct Ansatz, which is best constructed by exploiting the 
symmetries of our problem. Let us impose, first of all, our solution to be invariant 
under time-reversal t — > —t combined with L — > — L and R — > —R, under which also 
the CS term is invariant. This transformation reduces, in static configurations, to a 
sign change of the temporal component of L and R and of the spatial components 
of L and R. We can therefore consistently put them to zero. We also use parity 
invariance ({L <-» R, x <-» — x}) to restrict to configurations for which £j(x, z, t) — 
—Ri(—x,z,t), L 5 fi(x,z,t) — i?5,o(— x, z, t) and analogously for L, R. We impose, 
finally, invariance under "cylindrical" transformations [27], i.e. the simultaneous 
action of 3D space rotations x a o a — > 6^x a <r a 6, with 9 S SU(2), and vector SU(2) 
global transformations L, R — > 9 (L, R) #t. An equivalent way to state the invariance 
is that a 3D rotation with 6 acts on the solution exactly as an SU (2) vector one in 
the opposite direction (i.e. with 6') would do. The resulting Ansatz for the static 
solution (which we denote by "barred" fields) is entirely specified 4 real 2D fields 



1 



Rj (x, z) = Ai(r, z)x a Xj + -e ajk x k — e 

%(*,z) =A 2 {r, z)x a , 

— s(r, z) 
aR (x,z) = , 



"j 



(1.24) 



where r 2 = ^2 i x t x' 1 , x 1 = x t /r, e^ x ' v ^ is the antisymmetric tensor with e^ 1,2 ' = 1 
and the "doublet" tensors A^ 1 ' 2 ) are 



x a x b - S ab 



(1.25) 
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Substituting the Ansatz in the topological charge Eq. (|1.19[) we find 



B 



1 

2^ 



dr 



dp,(-i<f>*D p <p + h.c) + Afis 



(1.26) 



where 



= {r, z}, Ap = {Ai, A 2 }, A^ its field-strenght, <j) — 4>i + i(j} 2 and the 
covariant derivative will be defined in Eq. (|1.35p . The charge can be written, as it 
should, as an integral over the ID boundary of the 2D space. Finite-energy regular 
solutions with B = 1 which obey Eqs. (|1.3[) and (|1.4p must respect the following 
boundary conditions: 




zuv 




(1.27) 



and 



f <f*/r -> A 1 
0- + <h)/r- 
A 2 = 
s = 



A 2 = i 
s = 



(1.28) 



Solutions of the EOM with the required boundary conditions exist, and have been 
obtained numerically in Ref. [13] using the COMSOL package [28] (see Appendix 
for details). The 2D energy density of this solution is given in Fig. 11.11 



1.3.3. Zero-Mode Fluctuations 

Let us now consider time-dependent infinitesimal deformations of the static solu- 
tions. Among these, the zero-mode (i.e. zero frequency) fluctuations are partic- 
ularly important as they will describe single-baryon states. Zero-modes can be 
defined as directions in the field space in which uniform and slow motion is permit- 
ted by the classical dynamics and they are associated with the global symmetries 
of the problem, which are in our case U(2)v and 3-space rotations plus 3-space 
translations. The latter would describe baryons moving with uniform velocity and 
therefore can be ignored in the computation of static properties like the form factors. 
Of course, the global U(l)v acts trivially on all our fields and the global SU(2)y 
has the same effect as 3-space rotations on the static solution (|1.24|) because of 
the cylindrical symmetry. The space of static solutions which are of interest for us 
is therefore parametrized by 3 real coordinates -denoted as collective coordinates- 
which define an SU (2) matrix U. 

To construct zero-modes fluctuations we consider collective coordinates with 
general time dependence, i.e. we perform a global SU(2)v transformation on the 
static solution 



Rfi(x,z;U) = UR^z)U^ , i? (x,z;C/) = J2o(x,*), 



(1.29) 
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but we allow U = U(t) to depend on time. It is only for constant U that Eq. (|1.29|) 
is a solution of the time-dependent EOM. For infinitesimal but non-zero rotational 
velocity 



Eq. (jl.29|) becomes an infinitesimal deformation of the static solution. Along the 
zero-mode direction uniform and slow motion is classically allowed, for this reason 
our fluctuations should fulfill the time-dependent EOM at linear order in K provided 
that dK/dt = 0. 

From the action (|1.9| and (jl.lpp the following EOM are derived 



We only need to specify the EOM for one chirality since we are considering, as 
explained in the previous section, a parity invariant Ansatz. We would like to find 
solutions of Eq. (|1.30p for which Rp, and Rq are of the form p.290 : it is easy to see 
that the time-dependence of U in Eq. Q1.29[l acts as a source for the components 
Rq and Rp, which therefore cannot be put to zero as in the static case. Notice that 
the same happens in the case of the 4D skyrmion [4], in which the temporal and 
spatial components of the p and to mesons are turned on in the rotating skyrmion 
solution. Also, it can be shown that Eq. (|1.30p can be solved, to linear order in K 
and for dK/dt = 0, by the Ansatz in Eq. p.29[) if the fields i?o and Rp are chosen 
to be linear in K . Even though K must be constant for the EOM to be solved, 
it should be clear that this does not imply any constraint on the allowed form of 
the collective coordinate matrix U(t) in Eq. (| 1 . 29[) . which can have an arbitrary 
dependence on time. What we actually want to do here is to find an appropriate 
functional dependence of the fields on U(t) such that the time-dependent EOM 
would be solved if and only if the rotational velocity K = —iU^dU/dt was constant. 

In order to solve the time-dependent equations (|1.30|) we will consider a 2D 
Ansatz obtained by a generalization of the cylindrical symmetry of the static case. 
The Ansatz for Rp and Rq is specified by Eq. (|1.29|l in which the static fields are 
given by Eq. (jl.24|) . Due to the cylindrical symmetry of the static solution the fields 
in Eq. (|1.29p are invariant under 3D space rotations x a <j a — » 9'x a <J a 9 combined 
with vector SU{2) global transformations L,R — ► 9{L,R)9^ if U also transforms 
as U — ► 9^U9. We are therefore led to consider a generalized cylindrical symmetry 
under which k a also rotates as the space coordinates do. Compatibly with this 
symmetry and with the fact that Rq and Rp must be linear in K we write the 



K = k a a a /2 = -illUU/dt , 




(1.30) 



April 15, 2009 9:39 World Scientific Review Volume - 9.75in x 6.5in baryons 



Baryon physics in a five- dimensional model of hadrons 13 

Ansatz as 

i?o(x, z; U) = UR {x, z; K)U^ + i Ud U^ , # A (x, z; U) = z; K) , (1.31) 

where 

i?o(x, z- K) = x( x) (r, z)k b A {x ^ ab + v(r, z){k ■ x)x a 

a%(x, z- K) = f^y^- (k l - (k • + B t (r, z){k ■ x)x l + Q(r, z)e ibc k b x c 

ai?5(x, z; K) = B 2 {r, z){k ■ x) 

(1.32) 

It must be observed that our Ansatz has not fixed the 5D gauge freedom com- 
pletely; its form is indeed preserved by chiral SU(2)L t R gauge transformations of 
the form g R = U(t) ■ g ■ W(t) and g L = U(t) ■ g^ ■ U\t) with 

g = exp[ia(r, z)x a a a /{2r)] , (1.33) 

under which the 2D fields <f>M and X(x) defined respectively in Eq. Q1.24p and 
(|1.32p transform as charged complex scalars. The fields transform as gauge 
fields. There is also a second residual U(l) associated with chiral U(1)l,r 5D 
transformations of the form g^ = g and g^ = g' with 



g = exp 



i/3(r,z) 



(1.34) 



Under this second residual £7(1) only — {Bi, B 2 } and p transform non trivially; 
Bfi is a gauge field and p a Goldstone. In order to make manifest the residual 
gauge invariance of the observables we will compute we introduce gauge covariant 
derivatives for the </>, \ an( i P fields 

{Dp,cj>) {x) = df,(t) {x) + e ixy) Af,(/) {y) 

{D-uX)( x) =d fi x {x )+e {xv) A fiX(y) . (1.35) 
DfiP = dfip - Bfi 

At this point it is straightforward to find the zero-mode solution. The EOM 
for the 2D fields can be obtained by plugging the Ansatz in Eq. (jl.30[) . while the 
conditions at the IR and UV boundaries are derived from Eq. (|1.3[) and (|1.4[) . 
respectively. The boundary conditions at r = are obtained by imposing the 
regularity of the Ansatz, while those for r — > oo come from requiring the energy of 
the solution to be finite and B = 1. Also in this case, numerical solutions can be 
obtained with the methods discussed in the appendix. The reader not interested 
in detail can simply accept that a solution of Eq. (|1.30p exists and is given by our 
Ansatz for some particular functional form of the 2D fields which we are able to 
determine numerically. In the rest of the paper the 2D fields will always denote this 
numerical solution of the 2D equations. 
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1.3.4. The Lagrangian of Collective Coordinates 

The collective coordinate matrix U(t) will be associated with static baryons. The 
classical dynamics of the collective coordinates is obtained by plugging Eqs. (|1.29p 
and ljl.31jl in the 5D action. One finds S[U] = J dtL where 

(1.36) 



L 



-M + — k a k a . 
2 



The mass M and the moment of inertia A are given respectively by 



M = 8%M 5 
2 r 



and 



dr 



dz < a(z) 



\D n 



-r 2 A 2 - 
4 w 



^(i-i« 2 



d ii (-i<f>*D i? <f> + h.c) + Aj ls 



(1.37) 



1 



A = 16itM 5 - I dr 
3 J{ 

„2 



3 JO -I :r\ 



dz < a(z) 



- (DppY - r 2 (d- a QY - 2Q 2 - —BppBpB 



+r 2 (Djix) +ir(dav) + (x(x)X(x) + v 2 ) (l + (j)( x )(t)( x )) ~ 4v(j) (x) X( x ) 



- 7 Z 



- 2e^D n p X{x) (D p <t>) (x) + 2e^d p (rQ) X { x )^ y) 



9 (x)<P(x) 



1) +rQe* p An i? 



WD 



**> (£) 



(1.38) 

The numerical values of M and A are easily computed, once the numerical solution 
for the 2D fields is known. Using the best-fit values of the parameters we find 
M = 1132 MeV and 1/A = 227 MeV. 

Let us give some more detail on this theory. For now we proceed at the classical 
level and we will discuss the quantization in the next section. Our lagrangian can 
be rewritten as 



L = —M + ATr 



= -M 



(1.39) 



where we have parametrized the collective coordinates matrix U asU — ugl+i Ui<7 1 , 
with J2i u i 2 — 1- The lagrangian (|1.39[) is the one of the classical spherical rigid 
rotor. The variables {uQ,Ui\ are restricted to the unitary sphere S , which is 
conveniently parametrized by the coordinates q a = {x, (fiifo} -which run in the 
x S [—1, 1], 0i G [0, 2ir) and 4>2 € [0,27r) domains- as 



Ul +iu 2 = Zl = ^I_£ e <*, u + iu 3 = z 2 = ^j^^e^ 2 , (1.40) 

where we also introduced the two complex coordinates Z\£. We can now rewrite 
the Lagrangian as 



1 + x 



L 



-M + 2Xg a(3 q a q 13 , 



(1.41) 
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where g is the metric of S 3 which reads in our coordinates 



ds 2 = g aP dq a dqP = \—^—dx 2 + i-? dfa 2 + i±£ dfo 2 . (1.42) 



4 1 -x 2 



The conjugate momenta are p a — dL/dq a = A\g a pqP and therefore the classical 
Hamiltonian is 



H c = M + l.p a g<*0( q )pp. (1.43) 



It should be noted that the points U and — U in what we denoted as the 
space of collective coordinates actually describe the same field configuration (see 
Eq. Q1.29ll.31[l ). The SU(2) = S3 manifold we are considering is actually the uni- 
versal covering of the collective coordinate space which is given by S3 /Z% . This will 
be relevant when we will discuss the quantization. 



1.3.5. Skyrmion Quantization 

We should now quantize the classical theory described above, by replacing as usual 
the classical momenta p a with the differential operator —id/dq a acting on the 
wave functions /(<?). Given that the metric depends on q, however, there is an 
ambiguity in how to extract the quantum hamiltonian H q from the classical one 
in Eq. (|1.43[) . This ambiguity is resolved by requiring the quantum theory to have 
the same symmetries that the classical one had. At the classical level, we have an 
SO(A) ~ SU(2) x SU(2) symmetry under U -> U ■ 0t a nd U -> g ■ U with 6,g G 
SU(2). These correspond, respectively, to rotations in space and to isospin (i.e. 
global vector) transformations, as one can see from the Ansatz in Eqs. Q1.29ll.3ip . 
This is because K is invariant under left multiplication by g, and that the Ansatz 
is left unchanged by performing a rotation x a a a — > 0^x a a a 9 and simultaneously 
sending U — > U ■ 9. The spin and isospin operators must be given, in the quantum 
theory, by the generators of these transformations on the space of wave functions 
f(q) which are defined by 



[S a ,U] = Ucr a /(2), [I a ,U] = -a a /(2)U. (1.44) 
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After a straightforward calculation one finds 



S 3 
S+ 

s~ 
I 3 -. 
1+ 



"5< 
1 

1 

7T 



% x + 

i(0i+0 2 ) 
-i(0i+0 2 ) 



1 l-x 



1 + x 



<;VT 



1 1 + x, 



1 l-x 



~2 ^! ~ 

1 r i(<t>i-<h) 
V2 



I- = -=€ 

V2 



-i(<l>l-<j>2) 



iyl — x 2 d x - 
i\/ 1 — x 2 d x 



l-x 



1 1 + x, 



1 + x 



1 l-x 



(1.45) 



l-x 



1 + x 



1 l+x a 1 l-x a 



where the raising/lowering combinations are S f± = (S* 1 ± iS 2 )/y/2. 

The operators in Eq. (|1.45p should obey the Hermiticity conditions (S 13 ) = S 3 , 
{S + y = S~, and analogously for the isospin. In order for the Hermiticity conditions 
to hold we choose the scalar product to be 



(A\B) 



d 3 qJgfA ] {q)fB{q)- 



(1.46) 



where ^fg = 1/4 in our parametrization of S3. The reason why this choice of the 
scalar product gives the correct Hermiticity conditions is that S a and I a (where 
a = 1, 2, 3) can be written as X a d a with X a Killing vectors of the appropriate ^3 
isometries. The Killing equation V a X^ + V pX a = ensures the generators to be 
Hermitian with respect to the scalar product (|1 .46|) . 

Knowing that the scalar product must be given by Eq. (|1.46p greatly helps in 
guessing what the quantum Hamiltonian, which has to be Hermitian, should be. 
We can multiply and divide by ^fg the kinetic term of H c and move one yfg factor 
to the left of p a . Then we apply the quantization rules and find H 



(1.47) 



8A^5 vv ^ " 8A 
which is clearly Hermitian. We can immediately show that H q commutes with 
spin and isospin, so that the quantum theory is really symmetric as required: a 
straightforward calculation gives indeed 



M 



1 



-S z 



M 



1 



(1.48) 



2A 2A 

It would not be difficult to solve the eigenvalue problem for the Hamiltonian 
(|1.47p . but in order to find the nucleon wave functions it is enough to note that 
the versor of n-dimcnsional Euclidean space provides the n representation of the 

f The last equality holds because H q is supposed to be acting on the wave functions, which are 
scalar functions. 
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SO(n) isometry group. In our case, n = 4 = (2, 2), which is exactly the spin/isospin 
representation in which nucleons live. It is immediately seen that z\, as defined in 
Eq. p.40[) . has S 3 = I 3 = 1/2. Acting with the lowering operators we easily find 
the wave functions 



\p T> 
\p i) 



l 



-Zi . 



-Z2 ■ 



\n T> 
\n |) 



-zi 



(1.49) 



which are of course normalized with the scalar product (|1.46|) . The mass of the 
nucleons is therefore E = M + 3/ (8 A). 

Notice that the nucleon wave functions are odd under U — > —U, meaning that 
they are double- valued on the genuine collective coordinate space Sz/Zz- This 
corresponds, following [29], to quantize the skyrmion as a fermion and explains 
how we could get spin-1/2 states after a seemingly bosonic quantization without 
violating spin-statistic. 

Let us now summarize some useful identities which will be used in our calcula- 
tion. First of all, it is not hard to check that, after the quantization is performed 
the rotational velocity becomes 

1 



-i Tr 



U^Ua a 



A 



S a , 



and analogously 



iTr 



UU' i a a 



1- 



(1.50) 



(1.51) 



Other identities which we will use in our calculations are 
(Tr [Ua b tfa a ] = -^S b I a ) , 

(Tr [Ua b x b (k-x)tfa a ] = -^I a 



(1.52) 



where the VEV symbols (...) mean that those are not operatorial identities, but they 
only hold when the operators act on the subspace of nucleon states. Notice that the 
second equation in (II. 52ft is implied by the first one if one also uses the commutation 
relation (fHH]) . Eq. (fl~5T|) and the fact that, on nucleon states, ({S a ,S 1 } = S al /2). 



1.3.6. The Nucleon Form Factors 

The nucleon form factors parametrize the matrix element of the currents on two 
nucleon states. For the isoscalar and isovector currents we have 



<W/(p')l^(0)W(p)> =«/(?/) 



2M 



N 



(N f {p')\J» a (Q)\Ni{p))=u f {p') 



2M N 



(2i°)tti(p),(1.53) 
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where the currents are defined as J£ = J£ + J a L and J s = 1/3 (j R + J L ) in terms 
of the chiral ones. In the equation above q = p' — p is the 4-momentum transfer, 
Ni and Nf are the initial and final nucleon states and Ui(p), Uf(p') their wave 
functions, I a = a a /2 is the isospin generators and a^ v = i/2[7 M , 7"]. For the axial 
current J\ = Jjj — J£ we have 



(Nf(p')\J%JQ)\Ni(p)) =u f {p')G A (q*) 



2M N 

7/i —<f 

T 



^I a u f (p). (1.54) 



Exact axial and isospin symmetries, which hold in our model, have been assumed 
in the definitions above. 

In our non-relativistic model the current correlators will be computed in the 
Breit frame in which the initial nucleon has 3-momentum —q/2 and the final +<f/2 

(i.e. p» = {E,-q/2) and p'^ = (E,q/2), and q 2 = - f , with E = ^M^ + q 2 /A). 
Notice that the textbook definitions in Eqs. (|1.53ll.54p involve nucleon states which 
are normalized with y/2E; in order to match with our non-relativistic normalization 
we have to divide all correlators by 2M^. The vector currents become 

(JW/2)|jg(0)|JV«(-#2)) = G s E (f) X \xi, 

(N f (q/2)\Jm\Ni(-m) ^^^(Sxg)^, 
(N f (q/2)\Jv a (0)\Ni(-m) = G v E {f)x\{2I a ) X i> 
(N f tf/2)\Jv(Wi(-W2)) - *^^x}2(S ? xqT(2r) Xl , (1.55) 
where we defined 

G S E V (-1 2 ) = *f ' V) + T^rtf' V) . Gf} V (-q 2 ) = if ' V) + *f V) , 

and used the definition (S x q) % = e n S 3 q k . The nucleon spin/isospin vectors of 
state XiJ are normalized to x^X — !• F° r t ne axial current we find 

(^72)1^(0)1^(^72)} = x\^G A (f)2S T ^- X i, 
{N f {q/2)\jYmNr{-q/2)) = (1.57) 

where St = S — q S ■ q is the transverse component of the spin operator. 

It is straightforward to compute the matrix elements of the currents in posi- 
tion space on static nucleon states. Plugging the Ansatz 2411. 2911. 3211. 3"T]) in 
the definition of the currents (|1.13j) and performing the quantization one obtains 
quantum mechanical operators acting on the nucleons. The matrix elements are 
easily computed using the results of sect. 3.1. We finally obtain the form factors by 
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taking the Fourier transform and comparing with Eqs. (11.5511.5711 . We have@ 



E = _ 6^I I dr r Jo(qr) (a(z)d z s) 



uv 



3A 



drr 2 j (qr) a(z) (d z v - 2 (-D z x)( 2 )) 



uv 

3A 7 Q r 
v MnNq 2 jijqr) ( . w , \ 

where j' n are spherical Bessel functions which arise because of the Fourier transform. 



Ga = I <lr i 



a(z) JiM ( (jD ^ )(i) _ r A zr ) - a(z) j (qrj 



.58) 

(7V 



1.4. Properties of baryons: Results 

In this section we will present our results. After discussing some qualitative features, 
such as the large- N c scaling of the form factors and the divergences of the isovector 
radii due to exact chiral symmetry, we extrapolate to the physically relevant case 
of iV c = 3 and perform a quantitative comparison with the experimental data. 
Consistently with our working hypothesis that the 5D model really describes large- 
N c QCD we find a 30% relative discrepancy. 

Large-N c Scaling 

As explained in sect. 2.1, all the three parameters a, 7 and L of our 5D model should 
scale like N®, Eq. (|1.17p . in order for the large- N c scaling of meson couplings and 
masses to be correctly reproduced. This implies the following scaling for the baryon 
observables. First, we notice that the solitonic solution is independent of N c given 
that M5 factorizes out of the action and does not appear in the EOM. This implies 
that the radii of the soliton does not scale with N c , while the classical mass M and 
the moment of inertia A scale like N c . Using this we can read the A^-scaling of the 
electric and magnetic form factors from Eq. (|1.58l) : 

G S E ~N C , G V E ~N°, ^~iV c °, ^f-N c . (1.59) 

In large- N c QCD the baryon masses scale like A^ c [31], as in our model. The 
matrix elements of the currents on nucleon states are also expected to scale like 
N c , even though cancellations are possible [32]. The radii, therefore, must scale 

s It is quite intuitive that the form factors can be computed in this way. Given that solitons are 
infinitely heavy at small coupling, in the Breit frame they are almost static during the process of 
scattering with the current. To check this, however, we should perform the quantization of the 
collective coordinates associated with the center-of-mass motion, as it was done in [30] for the 
original 4D Skyrme model. 
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like AT? as we find and, looking at the definition (|1.55p . G S ^ V and /M^ should 
both scale like N c up to cancellations. It is very simple to understand why, both 
in QCD and in our model, there must be a cancellation in G\. Remembering 
that the temporal component of the current at zero momentum gives the conserved 
charge and looking at the definitions (|1.55p , one immediately obtains G\{fS) = 1/2 
because the skyrmion, as the nucleon, is in the 1/2 representation of isospin. This 
condition is respected by our model as it is implied by the EOM, and fulfilled to great 
accuracy (0.1%) by the numerical solution. Similarly we find at zero momentum 
G% (0) = N c /6 as required for a bound-state made of N c quarks of U(l)v charge 1/6 
each (in our conventions). Also this condition is implied by the EOM and verified 
by the numerical solution. 

Concerning the second cancellation, i.e. Gf f /Mjv ~ N®, we are not able to prove 
that it must take place in large- N c QCD as it does in our model. We can, however, 
check that it occurs in the naive quark model, or better in its generalization for 
arbitrary odd N c = 2 k + 1 [33] . In this non-relativistic model the Nucleon wave 
function is made of 2k + 1 quark states qi, 2k of which arc collected into k bilinear 
spin/isospin singlets while the last one has free indices which give to the Nucleon 
its spin/isospin quantum numbers. Of course, the wave function is symmetrized in 
flavor and spin given that the color indices are contracted with the antisymmetric 
tensor and the spatial wave function is assumed to be symmetric. The current 
operator is the sum of the currents for the 2k + 1 quarks, each of which will assume 
by symmetry the same form as in Eq. (| 1 . 55|) . If Si ,2 and I\ t 2 represent the spin and 
isospin operators on the quarks the operators Si + S2 and I\ + I2 will vanish on 
the singlet combination of the two quarks, but Si I\ + S2I2 will not. The k singlets 
will therefore only contribute to G§, G\j and Ga, which will have the naive scaling, 
while for the others we find cancellations. 

A detailed calculation can be found in [34] where, among other things, the 
proton and neutron magnetic moments and the axial coupling are computed in the 
naive quark model. The magnetic moments are related to the form factor at zero 
momentum as fiy/l^N = G\j{G) and hs/hn = Gf^O) where /j,n = 1/(2Mat) is the 
nuclear magneton and 2 fiy — fi p — 2 /is = fi p + /i„. In accordance with the 
previous discussion, the results in the naive quark model are 2^,$ = M« + an d 
2/iy = 2fc/3(/i M — fid), where \i u ,d are the quark magnetic moments, while for the 
axial coupling one finds §a — Ga(0) = 2fc/3 + 1 which scales like N c as expected: 

9A = T + i- (L60) 

Notice that for N c — 3 the subleading term in the l/A^-expansion represents a 60% 
correction. We have of course no reason to believe that such big corrections should 
persist in the true large- N c QCD; this remark simply suggests that "large" 1/N C 
corrections to the form factors are not excluded. 
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Divergences in the Chiral Limit 

It is well known that in QCD the isovector electric (r% v ) and magnetic (r 2 M v ) 
radii which are proportional, respectively, to the q 2 derivative of G\ and at zero 
momentum, diverge in the chiral limit [35]. In our model, as in the Skyrme model, 
divergences in the integrals of Eq. (|1.58p which define the form factors are due, as 
in QCD, to the massless pions. If all the fields were massive, indeed, any solution to 
the EOM would fall down exponentially at large r while in the present case power- 
like behaviors can appear. These power-like terms in the large-r expansion of the 
solution can be derived analytically by performing a Taylor expansion of the fields 
around infinity (1/r = 0), substituting into the EOM and solving order by order in 
1/r. The exponentially suppressed part of the solution will never contribute to the 
expansion. This procedure allows us to determine the asymptotic expansion of the 
solution completely, up to an integration constant j3. Substituting the expansion 
into the definitions of the form factors (|1.58|) one gets 



1 



oc (3 A J dr-jjoiqr) + 
G V E <x (3 2 J dr± 3o (qr) + 



G S M <x f3 3 



dr 



dr 



1 jijqr) 

r 5 qr 
1 jijqr) 
r 2 qr 



(1.61) 



All the form factors are finite for any q, including q = 0. The electric and magnetic 
radii, however, are defined as 



<4,M> 



G 



dG EM {q 2 ) 



GeM? 2 = °) 



dq"' 



(1.62) 



q 2 =0 



and taking a q 2 derivative of Eqs. (|1.6ip makes one more power of r 2 appear in the 
integral. It is easy to see that the scalar radii are finite, while the vector ones are 
divergent. For the axial form factor G a we find 



Ga oc I dr 



qr 



-Vt^ 3 )>"/''' 



(1.63) 



The integral in Eq. (|1.63[) is convergent for any q ^= but, however, it is not 
uniformcly convergent for q — * 0. The leading 1/r term in Eq. (|1.63|) is indeed 
given by I(q) — [3 J °° dr (1/r) (3ji(qr)/ (qr) — jo(qr)), which is independent of q 
and equal to /3/3, while the argument of the integral vanishes for q — > so that 
exchanging the limit and integral operations would give the wrong result 1(0) = 0. 
To restore uniform convergence and obtain an analytic formula for gA one can 
subtract the I(q) term from the expression in Eq. (|1.58p for Ga and replace it with 
/3/3. Rewriting the axial form factor in this way is also useful to establish that the 
axial radius, which seems divergent if looking at Eq. (j 1 . 63[) . is on the contrary finite. 
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The I(q) term, indeed, does not contribute to the q 2 derivative and the ones which 
are left in Eq. (|1.63[) give a finite contribution. 

We have found, compatibly with the QCD expectation, that all the form fac- 
tors and radii are finite but the isovector ones. Notice that the structure of the 
divergences is completely determined by the asymptotic large-r behaviour of the 
solution, and not by its detailed form (i.e., for instance, by the actual value of the 
integration constant (3 which depends on the entire solution). Our model coincides, 
in the IR, with the Skyrme model, therefore the asymptotic behaviour of the cur- 
rent densities is expected to be the same in the two cases. This explains why we 
obtained the same divergences as in the Skyrme model. 

Pion Form Factor and Goldberg er-Treiman relation 

It is of some interest to define and compute the pion-nucleon form factor which 
parametrizes the matrix element on Nucleon states of the pion field. In the Breit 
frame (for normalized nucleon states) it is 

(NftfMl^WlNii-fffl) = -^ jj ^G IiNv tf)x* f (2S i )q i {2I°)xi , (1.64) 

where 7r°(x) is the normalized and "canonical" pion field operator. The field is 
canonical in the sense that its quadratic effective lagrangian only contains the canon- 
ical kinetic term L2 = l/2((?7r Q ) 2 , or equivalently that its propagator is the canonical 
one, without a non-trivial form factor. With this definition, Gnntt is the vertex 
form factor of the meson-exchange model for nucleon-nucleon interactions [36] and 
corresponds to an interaction 

Cnn* = i {G NN „{n)ir a )N^ l5 (2I a )N . (1.65) 

On-shell, the form factor reduces to the pion-nucleon coupling constant, Gnntt{0) = 
gNN-K, whose experimental value is Qnntt = 13.5 ± 0.1. 

The pion field which matches the requirements above is given by the zero-mode 
of the KK decomposition. In the unitary gauge d z {a(z)A§) = 0, where Am = 
(Lm — i?Af)/2, and for AdS$ space, one has 

A^\x,z) = -±-^--K a (x)o- a , (1.66) 

where F^ is given in Eq. (|1 . 14|) . Gauge-transforming back to the gauge in which 
our numerical solution is provided and using the Ansatz in Eqs. (11.2411. 29p we find 
the pion field 

tt q = — ^ / dzA 2 (r, z)x b Ti [Ua b U^a a ] . (1.67) 



h Nucleon scattering, in our model, is a soliton scattering process and we have no reason to believe 
that it can be described by meson-exchange, i. e. that contact terms are suppressed. Therefore, we 
will not attempt any comparison of our form factor with the one used in meson-exchange models. 
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Table 1.2. Prediction of the nucleon observables with 
the microscopic parameters fixed by a fit on the mesonic 
observables. The deviation from the empirical data is 
computed using the expression (th — exp)/ min(|th|, |exp|), 
where th and exp denote, respectively, the prediction of 
our model and the experimental result. 





Experiment 


AdS 5 


Deviation 


M N 


940 MeV 


1130 MeV 


+20% 


MS 


0.44 


0.34 


-30% 




2.35 


1.79 


-31% 


9A 


1.25 


0.70 


-79% 




0.79 fm 


0.88 fm 


+11% 




0.93 fm 


oo 






0.82 fm 


0.92 fm 


+12% 




0.87 fm 


oo 






0.68 fm 


0.76 fm 


+12% 




-1.461 


-1.459 


+0.1% 



Taking the matrix element of the above expression and comparing with Eq. (|1.64p 
one obtains 

G NN ,{q 2 ) = -—M N F^q J drji{qr) J dzr 2 A 2 (r,z). (1.68) 

At q — ► the form factor Gnntt is completely determined by the large-r behavior 
of the field A 2 , given by A 2 — > f3/r 2 . We then find 

327T 

gNN-* = —MnF^L 2 . (1.69) 

By using Eqs. (33,34) of Ref. [2] which show that also qa is determined by the 
asymptotic behavior of the axial current, one finds 

gA = -^F 2 (3L 2 , (1.70) 

that, together with Eq. (|1.69|) . leads to the famous Goldberger-Treiman relation 
F-n9iTNN = M^gA- This relation, which is a consequence of having exact chiral 
symmetry, has been numerically verified to 0.01%. 

Comparison with Experiments 

Let us now compare our results with real- world QCD. We therefore fix the number 
of colors N c = 3 and choose our microscopic parameters to be those that gave the 
best fit to the mesonic quantities: 1/L ~ 343 MeV, M§L ~ 0.0165 and a ~ 0.94 
(7 ~ 1.23). The numerical results of our analysis and the deviation with respect 
to the experimental data are reported in tabic 11.21 Wc find a fair agreement with 
the experiments, a 36% total RMSE which is compatible with the expected size 
of 1/N C corrections. The axial charge gA is the one which shows the larger (80%) 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 '00 0.1 0.2 0.5 0.4 0.5 0.6 

<T [GeV 2 ) if (GeV-} 

Fig. 1.2. Scalar (left) and vector (right) electric form factors. We compare the results with the 
empirical dipole fit (dashed line) [4]. 




<? (GeV 3 ) if (GeV 2 ) 

Fig. 1.3. Normalized scalar (left) and vector (right) magnetic form factors. We compare the 
results with the empirical dipole fit (dashed line) [4]. 



deviation, and indeed removing this observable the RMSE decreases to 21%. We 
cannot exclude that, in a theory in which the naive expansion parameter is 1/3, 
enhanced 80% corrections to few observables might appear at the next-to-leading 
order. Nevertheless, we think that this result could be very sensitive to the pion 
mass and therefore could be substantially improved in 5D models that incorporate 
explicit chiral breaking. The reason for this is that gA is strongly sensitive to the 
large-r behavior of the solution (see the discussion following Eq. (|1.63[) ) which is in 
turn heavily affected by the presence of the pion mass. Notice that a larger value, 
gA — 0.99, is obtained in the "complete" model described in Ref. [4], a model with 
similar features to our 5D scenario and which includes a nonzero pion mass. This 
expectation, however, fails in the original Skyrme model, where the addition of the 
pion mass does not affect gA significantly [37] and one finds gA — 0.65. 

Table ITT21 also shows the proton-neutron magnetic moment ratio, /j, p /fj, n , which 
is in perfect agreement with the experimental value. This observable is the only one 
in the list that includes two orders of the 1/N C expansion. Indeed, due to the scaling 
/zy ~ N c and ^ s ~ N°, we have fip/[i n = -(/xy + fis)/(pv ~ Ms) - -1 - 2fJ.s/w- 

In figs. 11.21 11.31 and 11.41 we compare the normalized nucleon form factors at 
q 2 ^ with the dipole fit of the experimental data. The shape of the scalar and 
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0.0 0.2 0.4 0.6 0.S 1.0 0.0 0.2 0.4 0.6 0.8 L.O 

if (GeV) if (GeV 2 ) 

Fig. 1.4. Left: deviation of the ratio of proton and neutron magnetic form factors from the large 
N c value (solid line), compared with the dipole fit of the experimental data (dashed line). Right: 
normalized axial form factor (solid line) compared with the empirical dipole fit (dashed lines) [4] 
and with the experimental data taken from [38, 39] . 

axial form factors is of the dipole type, the discrepancy is mainly due to the error 
in the radii. The shape of vector form factors is of course not of the dipole type for 
small q 2 , but this is due to the divergence of the derivative at q 2 = 0. Including 
the pion mass will for sure improve the situation given that it will render finite the 
slope at zero momentum; it would be interesting to see if the dipole shape of these 
form factors is recovered in the presence of the pion mass. We also plot in the left 
panel of fig. 11.41 the deviation of ratio of the proton and neutron magnetic form 
factors from the large N c value which is given, due to the the different large- N c 
scaling of the isoscalar and isovector components, by G^, I (q)/G^ I (q) = — 1. Not 
only we find that this quantity is quite well predicted, with an error < 15%, but 
we also see that its shape, in agreement with observations, is nearly constant away 
from q 2 = 0. Also in this case corrections from the pion mass are expected to go in 
the right direction. 

1.5. Conclusions and outlook 

We have shown that five-dimensional models, used to describe meson properties of 
QCD, can also be considered to study baryon physics. Baryons appear in these 
theories as soliton of sizes of order 1 /m p stabilized by the presence of the CS term. 
We have reviewed the procedure to calculate the static properties of the nucleons 
that have shown to be in reasonable good agreement with the experimental data. 
This shows, once again, that 5D models provide an alternative and very promising 
tool to study properties of QCD in certain regimes. 

There are further issues that deserve to be analyzed. The most urgent one is the 
inclusion of a nonzero pion mass. As we have pointed out above, this will be crucial 
to calculate the isovector radii and, maybe, improve the prediction for gA- For this 
purpose we need to use a 5D model along the lines of Ref. [10, 11] where an explicit 
breaking of the chiral symmetry, corresponding to the quark masses, is introduced. 
We can also use this approach to study systems with high baryon densities, analyze 
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possible phase transitions or study the properties of nuclear matter. 
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A.l. Numerical Methods 

In this technical appendix we explain how the numerical determination of the soliton 
solution is performed. 



Equations of Motion and Boundary Conditions 

Let us first of all write down the EOM for the 2D fields which characterize our 
Ansatz in Eqs. (fOip . (TOTj) . (fT2T)|) and ([02]) . These can be obtained by plugging 
the Ansatz either directly in the 5D equations (|1.30[) or in the 5D action in Eqs. p.9[) 
and (|1.10[) . In the second case one gets the 2D action specified by Eqs. (|1.36[) . (|1.37[) 
and (|1.38[) and the EOM are obtained by performing the variation. In both cases 
one gets 

a(z) 



(o(«)D A < 



d p (r 2 a(z)A^) - a(z) (i<^D <j> + h.c.) + 1 Le^d- a (-) ( 
d- n (a(z)d»s) - [d fi {-i^D <f> + h.c.) + A^] = 



1) = , (A.l) 



for the fields which are already "turned on" in the static case. For the "new" fields 
which appear in the rotating skyrmion solution we have 

' d^{r 2 a{z)d fl v) - 2a(z) \v(l + |0| 2 ) - X ^ ~ fa* 

1 



+ ~fLe> J 



- ^Bpv + rQAjxp 



= 



D^{r 2 a(z)D- uX )+a{z) 2t^ - (1 + - jLe pp (D a 0) [id p (rQ) + D 9 p] = 



-d tl (r 2 a(z)d fl Q) - -a{z)Q 



l L -nu 



(W^iDtxV +h.c.) 



+\A nv {2v - X <P ~ fa f ) ~ ^D- aP d, (~) ] = Q ( A ' 2 ) 



(D^D.xV + h.c.) + -A-^x X ~ X0 f ) 



= 



r3 p (r 2 a{z)B^) + 2a(z)D p p 

+ 7 ie^{ [( X -v<t>)(D p <t>? 



h.c. 



(l-l^l 2 )^ 



2r 



'*(;)} 
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In order to solve numerically the EOM, they must be rewritten as a system of 
elliptic partial differential equations. This can be achieved by choosing a 2D Lorentz 
gauge condition for the residual Z7 (1) gauge fields 



0Mn = O, 



WBn = . 



(A.3) 



The equations for A 9 become J p = <9 P (r 2 aA^) = r 2 adfid^A" + d^ajA^ which 
is an elliptic equation and a similar result is obtained for Bp. 

The gauge condition needs only to be imposed at the boundaries, while in the 
bulk one can just solve the "gauge-fixed" EOM treating the two gauge field compo- 
nents as independent . The fact that the currents are conserved, dp J" = 0, implies 
indeed an elliptic equation for d^Ap which has a unique solution once the boundary 
conditions are specified. If imposed on the boundary, therefore, the gauge conditions 
are maintained also in the bulk. 

The IR and UV boundary conditions on the 2D fields follow from Eq. (II. 3p and 
Eq. (|1.4|) and from the gauge choice in Eq. (|A.3j) . They are given explicitly by 



z — Ziyi '■ 



and 



z = zijy 



0i =0 
d 2 (f>2 = 
Ai=0 
d 2 A 2 = 
I d 2 s = 



01 =0 

02 = -1 

Ax = 
d 2 A 2 = 
I s = 





(A.4) 



Xi =0 

X2 = -1 

v = -1 

10 = 




(A.5) 



The boundary conditions at r — oo have to ensure that the energy of the solution 
is finite; this means that the fields should approach a pure-gauge configuration. At 
the same time one has to require that the solution is non-trivial and its topological 
charge (Eq. (|l.I9p ) is equal to one. We have 




X7TZ / L 




z/L 




(A.6) 



The r — boundary of our domain requires an ad hoc treatment, given that the 
EOM become singular there. Of course this boundary is not a true boundary of 
our 5D space, but it represents some internal points. Thus we must require the 2D 
solution to give rise to regular 5D vector fields at r — and we must also require 
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the gauge choice to be fulfilled. These conditions are 




p/r -> B x 
SiBi = 
B 2 = 



(A.7) 



A. 1.1. COMSOL Implementation 

To obtain the numerical solution of the EOM we used the COMSOL 3.4 package [28] , 
which permits to solve a generic system of differential elliptic equations by the finite 
elements method. A nice feature of this software is that it allows us to extend the 
domain up to boundaries where the EOM are singular (i.e. the r = line), because 
it does not use the bulk equations on the boundaries, but, instead, it imposes the 
boundary conditions. 

In order to improve the convergence of the program and the numerical accuracy, 
one is forced to perform a coordinate and a field redefinition. The former is needed 
to include the r = oo boundary in the domain in which the numerical solution is 
computed. The advantage of this procedure is the fact that in this way one can 
correctly enforce the right behaviour of the fields at infinity by imposing the r = oo 
boundary conditions. A convenient coordinate change is given by 



where x is the new coordinate used in the program and c is an arbitrary constant. 
The domain in the x direction is now reduced to the interval [0, ctt/2]. The param- 
eter c has been introduced to improve the numerical convergence of the solution. A 
good choice for c is c ~ 10, which allows to have a reasonable domain for x and, at 
the same time, does not compress the solution towards x = 0. 

A field redefinition is needed to impose the regularity conditions at r = 
(Eq. (|A.7j) ). For this purpose we use the rescaled fields 



With these redefinitions, in the new coordinates, the r = boundary conditions 
read as 



In order to ensure the convergence of the program another modification is 
needed. As already discussed, to obtain a soliton solution with non- vanishing topo- 
logical charge we have to impose non-trivial boundary conditions for the 2D fields 




(A.8) 




(A.9) 




(A.10) 
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at r = oo (Eq. (|A.6[) ). It turns out that if imposing such conditions the program 
is not able to reach a regular solution. This is so because the r = oo boundary 
is singular and imposing non-trivial (though gauge-equivalent to the trivial ones) 
boundary conditions at a singular point spoils the regularity of the numerical solu- 
tion; the same would happen if the topological twist was located at r = 0. To fix 
this problem we have to perform a gauge transformation which reduces the r = oo 
conditions to trivial ones and preserves the ones at r = at the cost of introducing 
a "twist" on the UV boundary. For this, we use a transformation of the residual 
Z7(l) chiral gauge symmetry associated to SU(2)l : r (Eq. (|1.33[l ) with 



a(r,z) = (l-z/L)f(r), 
where /(r) can be an arbitrary function which respects the conditions 



l/(°o)-' 



and 



f" (0) = 
/"(oo)->0 ' 

is f(r) 



(A.ll) 



(A.12) 



-sin/(r) 
cos/(r) 
-1 



For c ~ 10 it turns out that a good choice for f(r 
gauge-fixing condition for Ap_ is now modified as 

d r A 1 +d z A 2 -(l-z/L)f"(r)=0 : 

the UV boundary conditions are given by 

xipi = sin/(r) 
(-1 + ^a) = -cos/(r) 
M = f(r) 
d z A 2 = 
s = 



and the r = oo constraints are now trivial 

(-l + a^a) = 1 
x : < d x Ai = 

^2=0 

s = 



2arctanr. The 



(A.13) 










(A.14) 



(A.15) 



whereas the r = and the IR boundary conditions are left unchanged. Notice that 
in the new gauge the EOM for A^ are modified in accord to Eq. (|A.13[) . however 
they are still in the form of elliptic equations. 
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